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We study an extended SU(N) single-impurity Kondo model in which the impurity spin is described 
by a combination of Abrikosov fermions and Schwinger bosons. Our aim is to describe both the 
quasiparticle-like excitations and the locally critical modes observed in various physical situations, 
including non-Fermi liquid behavior in heavy fermion systems in the vicinity of a quantum critical 
point. We identify the strong coupling fixed point of the model and study its stability within second 
order perturbation theory. Already in the single channel case and in contrast with either the pure 
bosonic or the pure fermionic case, the strong coupling fixed point is unstable against the conduction 
electron kinetic term as soon as the amount of Abrikosov fermions reaches a critical value. In the 
stability region, the partially screened, dressed impurity at site repels the conduction electrons on 
adjacent sites. In the instability region, the impurity tends to attract (TV — 1) conduction electrons 
to the neighboring sites, giving rise to a two-stage Kondo effect with additional screening. This 
result opens the route to the existence of an intermediate coupling fixed point, characterized by 
non-Fermi liquid behavior. 



I. INTRODUCTION 

Recent experiments in Heavy-Fermion compounds 
have shown the existence of a quantum phase transition 
from a magnetically disordered to a long-range magnetic 
ordered phase, driven by change in chemical composi- 
tion, pressure or magnetic fieldcl. For an extensive sur- 
vey of the experimental situatioa we refer the reader 
to the review article of Stewartcl. In a very unusual 
way, the behavior of the system in the disordered phase 
close to the quantum critical point (QCP) diffets-from 
that of a Fermi-liquid. For example CeCue- x Auja^u and 
(Cei^ x La x )Ru2Si^ present an antiferromagnetic tran- 
sition, respectively, at x c =0.1 and at x c = 0.08. While 
far from the QCP, the magnetically disordered phase is 
a Fermi liquid with a large effective mass, the tempera- 
ture dependence of the physical quantities in the disor- 
dered phase in the vicinity of the QCP is of non-Fermi- 
liquid like type. Typically, in CeCu$. gAuQ ia, the spe- 
cific heat C depends on T as C/T ln(T/T ), the 

magnetic susceptibility as \ ~ 1 — oty/T, and the T- 
dependent part of the resistivity as Ap ~ T instead 
of C/T ~ x ~ Const and Ap ~ T 2 as in the Fermi- 
liquid state. Once a long-range magnetic order is set 
up, the effect of a pressure or of a magnetic field is to 
drive the system back to a magnetically disordered phase 
with a Fermi-liquid behavior. The same type of behavior 
has been observed in other systems such as YbRh2Si2-, 
CeNi 2 GeM, CeCn 2 {Si 1 - x Ge x ) 2 B, Celn 3 , CePd 2 Si 2 B, 
and Ui- x Y x Pdji3. The associated breakdown of the 
Fermi-liquid theory poses fundamental questions about 
the possible formation of novel electronic states of matter 
with new types of elementary excitations resulting from 
the presence of strong correlations among electrons. 

On the theoretical side, two scenarios are in compe- 



tition to describe quantum phase transitions: either the 
itinerant magnetism scenario (i), or more recently pro- 
posed, the locally critical picture (ii). 

In the former case, (i), the quasiparticles still exist at 
the QCP and the theory focuses on the study of the low- 
lying, large- wavelength (low-w, low-g) fluctuations of the 
order parameter close to the transition. The calculations 
have b co j ji pe rformed within the renormalization group 
schemeOliBliia or in the self-consistent spiiw-fiuctuation 
theorytJ, and have been recently extendedEJ to the mi- 
croscopic model which is believed to describe the Heavy 
Fermions, the Kondo lattice. In all the cases, they lead 
to a $ -theory with an effective dimension d e ff = d + z 
where d is the spatial dimension and z is the dynamic 
exponent. In the experimental situations, d e ff is above 
its upper critical value equal to 4, since d is equal to 3 or 
2, and z varies from 2 to 3 depending whether the spin 
fluctuations are staggered or uniform. Hence the sys- 
tem is described by a Gaussian fixed point with anoma- 
lous temperature dependence of C/T and a — Ap/T but 
with predictions which cannot account for the non- Fermi- 
liquid behavior observed experimentally. 

The second scenario, (ii), has been motivated by the 
recent results obtained by inelastic neutron scattering ex- 
periments performed on CeCu^gAuQ^. The dynamical 
spin susceptibility x"(q, w) near the magnetic instabil- 
ity wavevector Q has-, been found to obey an anoma- 
lous uj/T scaling Iawll3 : t3 as a function of temperature: 
x"(Q,uj) ~ T~ a g(to/T) with an exponent a of order 
0.75. Moreover, such a u> and T dependence appear to 
stand over the entire Brillouin zone revealing in the 
bulk susceptibility too. This fact strongly suggests that 
the spin dynamics are critical not only at large length 
scalesE3E2l but also at atomic length scales contrary to 
what happens in the traditional itinerant magnetism pic- 
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ture, (i). From these results, one can deduce that local 
critical modes coexist with large-wavelength fluctuations 
of the order parameter implying a non-Gaussian fixed 
point beyond the <£> 4 -theory. Alternative theories to the 
spin-fluctuation scheme are needed to describe the lo- 
cal feature of the quantum critical point characterized 
by the simultaneous disappearance of the quasiparticles 
and the formation of local moments. In this direction, 
we will mention recent, calculations based on a dynami- 
cal mean field theoryEJEj (DMFT) which seem to lead to 
encouraging results concerning the scaling law variation 

of x"(q,w). 

A recent approach which has been developed in order 
to desimj^JJicJocal QCP, is based on a supersymmetric 
theoryEjSo'LJ in which the spin is described in a mixed 
fcrmionic-bosonic representation. The interest of the su- 
persymmetric approach is to describe the quasiparticles 
and the local moments on an equal footing through the 
fermionic and the bosonic part of the spin, respectively. 
It appears to be specially well-indicated in the case of 
the locally critical scenario in which the magnetic tem- 
perature scale T/v, and the Fermi scale Tk (the Kondo 
temperature) below which the quasiparticles die, vanish 
at the same point, 5c- 

An important aspect in the discussion of the break- 
down of the Fermi-liquid theory is related to the question 
of the stability of the strong coupling (SC) fixed point. 
Whereas all the issues presented previously concerning 
Heavy fermion systems have to do with properties of the 
lattice, the instability of the SC fixed point can be re- 
garded already by studying the single impurity problem. 

The traditional source of instability in the single im- 
purity Kondo model is the presence of several channels 
for the conduction electrons with the existence of two 
regimes, underscreened and the overscreened, with very 
different behaviors as we are about to recall. Indeed we 
will see that this is not the only possible source of insta- 
bility of the strong coupling fixed point. Recent works 
have shown that more general Kondo impurities of sym- 
metry group SU(N) may also lead to an instability of 
the SC fixed point already with one channel of conduc- 
tion electrons. 

In order to fix ideas, let us start with the antiferro- 
magnetic single-channel Kondo impurity model. It is 
well kn own . that within a renormalization group (RG) 
analysisElffiU'Ej, the flow takes the Kondo coupling J all 
the way to strong coupling. The weak-coupling beta 
function follows the renormalization group equation 



0(g) = 



dg(A) 
dA 



-9 



(1) 



where g = poJ and po is the density of states of con- 
duction electrons. The system flows a to strong coupling 
fixed point which is stable and the associated behavior 
of the system is that of a local Fermi liquid. 

The situation is rather different when one considers 
several channels for the conduction electrons. In the case 



of a spin 5 of symmetry group SU(2) in Kondo interac- 
tion with conduction electrons belonging to K different 
channels, Blandin and NozieresEa have shown that the 
multichannel Kondo model can lead to two very differ- 
ent situations depending how K compares to 25. Their 
calculation corresponds to a second order perturbation 
theory in the hopping amplitude, t, of the conduction 
electrons, around the strong coupling fixed point. They 
analyze their results by deriving an effective coupling, 
J e ff between the spin of the composite formed by the im- 
purity dressed by the conduction electrons in the strong 
coupling limit, and the spin of the conduction electron 
on the neighboring sites. They are then able to apply 
the same RG analysis to J e ff as indicated in Eq. (|). 
In the underscreened regime, when K < 25, the effec- 
tive coupling is found to be ferromagnetic and the strong 
coupling fixed point is stable. In the overscreened regime 
when K > 2S, the effective coupling is found to be anti- 
ferromagnetic and hence the strong coupling fixed point 
is unstable. The former K < 2S regime corresponds to 
the one-stage Kondo effect with the formation of an ef- 
fective spin, (5—1/2), resulting from the screening of the 
impurity spin by the conduction electrons located on the 
same site. The system described by the strong coupling 
fixed point, behaves as a local Fermi liquid. The insta- 
bility of the strong coupling fixed point obtained in the 
latter regime, K > 25, is associated with a multi-stage 
Kondo effect in which successively the impurity spin is 
screened by conduction electrons on the same site, and 
then the resulting dressed impurity is screened again by 
conduction electrons on the neighboring site and so forth. 
The instability of the strong coupling fixed point in the 
oversreened regime is an indication of the existence of an 
intermediate coupling fixed point which has been then 
investigatedEEfEdH by means of other methods. As it 
is well established now, the intermediate coupling fixed 
point leads to non- Fermi-liquid excitation spectrum with 
an anomalous residual entropy at zero temperature. 

It has recently put forward that other sources of in- 
stability of the SC fixed point may exist else but the 
multiplicity of the conduction electron channels. Recent 
works have shown that the presence of a more general 
Kondo impurity where the spin symmetry is extended 
from SU(2) to SU(N), and the representation is given 
by a L-shaped Young tableap, may also lead to an insta- 
bility of the SC fixed pointE-aEa already in the one-channel 
case. In the large N limit, Coleman et al. have found 
that the SC fixed point becomes unstable as soon as q 
(the number of boxes in the Young tableau along the 
first column), is larger than A^/2 whatever the value of 
25 (the number of boxes in the Young tableau along the 
first row) is. This result opens the route to the existence 
of an intermediate coupling fixed point with presumably 
non- Fermi-liquid excitation spectrum. The consideration 
of a L-shaped Kondo impurity fits in with the supersym- 
metry approach that we have evoked before since both 
spin operators and states can be expressed in terms of 
bosons and fermions. 
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At that point, it is worth noting that the supersym- 
metry theory, or specifically taking into consideration 
more general L-shaped Kondo impurities appears to offer 
valuable insights into the two issues raised by the break- 
down of the Fermi liquid theory that we have summarized 
above, i.e. both the existence of locally critical modes 
and the question of the instability of the SC fixed point. 
Somehow it seems that the choice of L-shaped Kondo im- 
purities captures the physics present in real systems, with 
the coexistence of the screening of the spin by conduction 
electrons responsible for the formation of quasiparticles, 
and the formation of a localized magnetic moment that 
persists and eventually leads to a phase transition as the 
coupling to other impurities becomes dominant. In the 
same way as large N expansions may provide insights 
into real systems even at finite value of the degeneracy, 
the study of more general impurities may enlighten the 
understanding of experimental situations. 

The aim of the paper is to study the extended 
SU(N) L-shaped single-impurity Kondo model in the 
one-channel case. We want to understand how the system 
behaves, not only as a function of the impurity param- 
eters, (25, q), but also as a function of the number of 
electrons, rid, available on neighboring sites, that is to 
say, of the filling. We find that as long as the bosonic 
component of spin is of order TV, there is a transition 
around the point where the fermionic component of the 
impurity is q = N/2. At this particular point, the energy 
shift is, to lowest order in perturbation theory around the 
strong-coupling fixed point, equal to (— 2t 2 / J), indepen- 
dently of the impurity parameters, q, S and N. When 
q < N/2, the low-energy physics corresponds to a stable 
strong-coupling fixed point. For q > N/2 the strong- 
coupling fixed point is unstable and anomalous behavior 
is expected, in particular, a two-stage quenching effect, 
as predicted by Ref. ^| This phase diagram is not acci- 
dental, but is due to the relation of the effective dressed 
impurity in the strong-coupling regime to the conduction 
electrons in neighboring sites, as our study of the depen- 
dence of the energy shifts on rid reveals. If q < N/2, 
the energy is minimized when the dressed impurity re- 
pels the electrons on the next site. That is, when = 1. 
At q = N/2, the energy shift is independent of rid- Fi- 
nally, the lowest energy shift for q > N/2 corresponds 
to a maximal rid, indicating the accumulation of conduc- 
tion electrons on neighboring sites, leading to a two-stage 
Kondo quenching. 

The rest of the paper is organized as follows. In Sec- 
tion [n], we introduce the model and the main features 
of the strong-coupling limit, where the electron kinetic 
term is neglected. In this limit the model is reduced to 
a single site problem, where the impurity is coupled to 
n c conduction electrons. We identify the ground state 
and the energies of the excited states with one more or 
one less conduction electron, which will play a role in the 
lowest order in perturbation theory. In section III, we 



discuss the stability of the strong coupling fixed point as 



resulting from the sign of the effective coupling J e ff be- 
tween the spin of the strong coupling composite and the 
spin of the conduction electrons on the neighboring site. 
The calculation is based on a second order perturbation 
theory in t and is performed for an arbitrary number rid 
of conduction electrons on the neighboring site. Section 
IV contains the discussion of the results. In the large 
N limit, we show how J e ff is derived from the energy 
shift difference between the symmetric and the antisym- 
metric configurations, and how the analysis of the rid 
dependence of the energy shift provides information on 
the nature of Kondo screening, with the realization of a 
two-stage Kondo effect when the SC fixed point becomes 
unstable. When the behavior of the system is controlled 
by the strong coupling fixed point, i.e. when q < N/2, 
the impurity in the ground state tends to repel electrons 
on neighboring sites. Once q > N/2, the repulsion be- 
comes attraction. We show how this feature is already 
present in the purely fermionic case, and is a consequence 
of the particle-hole symmetry (more details are contained 
in Appendix ^) . The fact that there is extra degeneracy 
in the supcrsymmetric impurity, due to the bosonic con- 
tribution, leads to the instability of the strong coupling 
fixed point as soon as q > N/2. We finish the section with 
a short discussion on the behavior of physical quantities 
in the different regimes. 

The appendices contain the technical details of the cal- 
culations. In Appendix |X| we outline the construction of 
three particle states with SU (3) symmetry, as an intro- 
duction to the group theoretical formalism used. Explicit 
expressions for the impurity states and the eigenstates 
of the model in the strong coupling limit are derived in 
Appendix ||. We also include a general presentation of 
the different representations of the spin, either bosonic, 
fermionic or L-shaped, as considered in the paper. We 
will show how in the latter case the spin operators and the 
impurity states are expressed in terms of fermion and bo- 
son creation and annihilation operators within two con- 
straints. Appendix |c| contains a calculation of the energy 
shift to the strong coupling fixed point to lowest order 
in perturbation theory, for the completely antisymmetric 
impurity. Since the ground state is a singlet, there is no 
splitting of levels. Nevertheless, the behavior of the en- 
ergy with the filling, rid, on the neighboring site shares 
many common features with the problem that we have 
studied. Finally, we include the details of the calcula- 
tion of the matrix elements needed in the second order 
perturbation theory calculation in Appendix 



II. THE MODEL AND ITS 
STRONG-COUPLING LIMIT 

A. SU(N) single-impurity Kondo model 

We consider a generalized, single-impurity, Kondo 
model with one channel of conduction electrons and a 
spin symmetry group extended from SU(2) to SU(N). 
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An impurity spin, S, is placed at the origin (site 0). In 
this article we will deal with impurities that can be real- 
ized by a combination of bosonic and fermionic operators, 
and are thus described by a L-s hap e d . representation in 
the language of Young tableaujSESEa, as illustrated in 
Fig. [I] (for details, see Appendix |b|). 



R = q — 



FIG. 1: Young Tableau description of an impurity with 
mixed symmetry, [2S, realized by a combination of 

fermions and bosons 

If 2S and q are the numbers of boxes along the first 
row and the first column, respectively, the representation 
is denoted by [2S, I 9-1 ]. Its degeneracyES is reported 
in Table 1. The conduction electrons transform under 
the fundamental representation of SU(N) and can be 
represented by Young tableaux made out of single boxes. 
The dimension of the fundamental representation is N 
which just means that each electron can be in one of N 
states of spin. 

The hamiltonian describing the model is written as 

H = £ £ k 4 ia c k , Q + J ]T S A £ c t (0)r^ C/3 (0) , (2) 



a,0 



where c^. is the creation operator of a conduction elec- 
tron with momentum k, and SU(N) spin index a — 



a,b,...,rN, 4,(0) 



is the creation oper- 



ator of a conduction electron at the origin, N$ is the 
number of sites, and T A g (A — 1, . . . , iV 2 — 1) are the 
generators of the SU (N) group in the fundamental rep- 
resentation, with Tt[t a t b ] = 5 AB /2. In the SU(2) 
case, t a = cr A /2, where {cr A } are the Pauli matri- 
ces. The conduction electrons interact with the impu- 
rity spin S A (A = 1, . . . , iV 2 — 1), placed at the origin, 
via Kondo coupling, J > 0. When the impurity is in 
the fundamental representation, we recover the Coqblin- 
Schrieffer modclExj describing conduction electrons in 
interaction with an impurity spin of angular momentum 
j, (N = 2j + 1), resulting of the combined spin and orbit 
exchange scattering. In our notation, a = j, b = j — 1,..., 

nv = -j- 



B. Strong-Coupling fixed point 

In the strong-coupling limit, the hamiltonian reduces 
to the local Kondo interaction term at site 

H = jJ2 sA Y, c «(°)T A cp(O), (3) 

A a, 13 



The ground state, \GS), is formed by binding the right 
amount of conduction electrons to the impurity in order 
to minimize the Kondo energy. Let us denote by Y (Fig. 
||) the representation of the n c conduction electrons cou- 
pled to the impurity, R that of the free impurity (Fig. ^), 
and Rsc the representation of one of the strong-coupling 
states resulting of the direct product R ® Y (cf. Fig. ||) 
(see Appendix |b] for details) . 



Y 



FIG. 2: Young tableau description of n c conduction electrons, 
localized at the impurity site 



2S 



N < —I 

c^ 

c 



Rsc 



FIG. 3: Young tableau description of the formation of the 
strong coupling ground state. We denote the presence of con- 
duction electrons at site by c. Notice that the first column 
in the Young tableau for Rsc is a singlet and can be removed. 

When N = 2, the Kondo energy can be written in 
terms of conserved quantities 

JS.£ C t(0)r Q/3C/3 (0) \GS} = 

a, (3 

J [ S sc {s sc + s R {s n + s y (s y + 1} ] | G5) ; 

where S(S + 1) is the eigenvalue of the Casimir operator 
S 2 for N = 2. The generalization to SU(N) is given by 

A a, 13 

~ [c2(Rsc) - C 2 {R) - C 2 (Y)] \GS) , (4) 

where C 2 {R) is the quadratic Casimir operator of the rep- 
resentation R, defined earlier, which commutes with all 
the generators of the group. For a representation given 
by a Young Tableau with rrij boxes in the j-th row until 
the row j = h, the eigenvalue C2({mj}) of the quadratic 
Casimir operator is 



C 2 ({m j }) 



Q(N 2 - Q) 
N 



TTij (rrij + 1 — 2j) 

i=i 
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where Q = J2j=i m j 1S the total number of boxesEa. Ta- 
ble H summarizes the expressions of the Casimir eigenval- 
ues for the impurities described in this work and for the 
conduction electrons, as well as the dimension of their 
spin representations. 

Minimization of the energy Eq. (^) , leads to a ground 
state with n c — (N — q) conduction electrons coupled to 
the L-shaped Kondo impurity ensuring partial screening. 
The resulting composite at site 0, with energy E , is made 
out of the impurity dressed by the conduction electrons 
in order to form a singlet along the first column. The 
associated Young tableau in the strong coupling regime 
is given in Fig. ^. Note that the first column of length 
N can be removed without changing the representation 
since it is a singlet. When the strong-coupling fixed point 
is stable, this corresponds to a one-stage Kondo effect in 
which the impurity is screened by the conduction elec- 
trons to form a bosonic (S — 1/2) impurity. 



Ground state 



quires an additional conduction electron at the impurity 
site, \GS+1), or it loses one conduction electron, \GS— 1). 
In the former case, \GS + 1), the spin of the additional 
conduction electron can be either symmetrically or an- 
tisymmetrically correlated with the spin of the impurity 
as schematized in Fig. [|. 



\GS + 1) 1 



\GS + 1) A = — 



Let us now write the expression of the fundamental 
state associated with this strong-coupling fixed point. 
The ground state is degenerate. The states in the multi- 
plet transform as a completely symmetric representation 
of SU(N), described by a Young tableau with (25 — 1) 
boxes, denoted by [25 — 1], (Fig. ||). We choose a realiza- 
tion of the impurity in terms of 25 bosonic operators and 
(q — 1) fcrmionic operators, which happens to be more 
convenient. We could have constructed impurity states 
with the same SU(N) symmetry using (25 — 1) bosons 
and q fermions (see Appendix |B|). We would like to em- 
phasize that all the results that we establish in this paper 
are independent of the operator representation which we 
choose to work with. The highest weight state is then 
written as 

1 M) 2S ~ X \ A ) . (5) 



D! 



with 



1 l q iN 

\A)^-A(b{(i[fi)( n 



4))io> 



(6) 



= iq+l 



7 ^ ^/ (2S + N - l)C q N _\ . 

Here, |A) transforms itself as a SU(N) singlet and it 
will be annihilated by any of the raising and lowering 
operators, T^A) = U ± \A) = ■■■ = 0. This "state" 
would describe the strong coupling ground state for a 
purely fermionic impurity. 



D. Excited States 

There are two types of excited states in the strong- 
coupling regime. Either the degenerate ground state ac- 



FIG. 4: Excited states, \GS + 1} S , \GS + 1} A , with an addi- 
tional conduction electron, n c = (TV — q + 1), respectively in 
the symmetric and antisymmetric configurations. 

In the limiting case of SU(2) spin, these two configura- 
tions correspond to a spin of the conduction electron that 
is either parallel or antiparallel to the impurity spin. In 
the general SU (TV) case, we will keep on speaking of sym- 
metric and antisymmetric configurations respectively. 

States with one less electron will be denoted by \GS— 1) 
and are represented by the Young tableau in Fig. ||. Let 



\GS-1) = (N-l) 



FIG. 5: Excited state, \GS 
electron n c = (TV — q — 1). 



1), with one less conduction 



us denote by A£f = Ef - E , AE^ = - E and 
AE2 = E2 — E0 the energy differences, with respect to the 
ground state energy, associated with these three excited 
states |G5 + 1) 5 , \GS+1) A and |G5-1>. Using the same 
casimirology method as presented at the beginning of this 
section for the determination of the ground state energy, 
we have summarized our results in Table [l| respectively 
for arbitrary N and in the large-N limit with (25 + q — 
1)/TV finite. 

One can check that the results in Table |n| coincide with 
(Eqs. [25-26]) in ref. |25|, within a TV/2 factor stemming 
from a different definition of the Kondo coupling J (cf. 
Eq.(l) of ref. ||) and of the Casimir (cf. Eq.(17) of ref. 
|25| ), and a change in the notations n*j — q and ni, = 25. 
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TABLE I: Dimension, d, and eigenvalues of the casimir operator, C2 for the symmetric, antisymmetic, L-shaped, and funda- 
mental representations studied in this paper. In the L-shaped case, Q = (2S + q—l) is the total number of boxes in the Young 
tableau, and Y' = (q — 2S) measures the row-column asymmetry 





Symmetric 
[25] 


Antisymmetric 
[1"] 


L-shaped 

[2S, I"' 1 ] 


Fundamental 
[1] 


d 


n 2S 

^N+2S-1 




( 2S \ r<2S ni- 1 

\2S+q-l j L 'N+2S-1^N-1 


N 


c 2 


^(2S(2S + N)(N-1)) 


±( q (N-q)(N + l)) 


%{N-Y' - Q/N) 
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TABLE II: Strong coupling excitation energies, AEf , AEi, 
and AE2, in the case of an L-shaped impurity, measured with 
respect to the ground state, of the states with one more con- 
duction electrons on site 0, coupled symmetrically and anti- 
symmetrically, respectively, to the dressed impurity on site 0, 
and of the state with one less electron. 





AE? 


AEJ 1 


AE 2 


Arbitrary N 


i(2S + N - q - Q/N) 


i(N-q-Q/N) 


Hq + Q/N) 


Large N limit 








(Q/N finite) 


i(N-q+2S) 


i(N-q) 


ii 



III. STABILITY OF THE STRONG COUPLING 
FIXED POINT 



We have identified the strong coupling fixed point in 
the previous section. For J — * 00, the lowest energy state 
corresponds to n c — (N — q) electrons partially screening 
the impurity at the origin, and free electrons in the other 
sites, unable to hop to the impurity site. 

In order to better understand the low-energy physics 
of the system, we should consider the finite Kondo cou- 
pling, allowing virtual hopping from and to the impu- 
rity site. These processes generate interactions between 
the composite at site and the conduction electrons on 
neighboring sites, that can be treated as perturbations 
of the strong coupling fixed point. Applying an analysis 
similar to that of Nozieres and BlandincJ to the nature of 
the excitations, we can argue whether or not the strong 
coupling fixed point remains stable once virtual hopping 
is allowed. 

We consider a system with an additional site next to 
the dressed impurity, filled with n& electrons. The ground 
state consists of two multiplets, with different symmetry 
properties. Once the hopping is turned on, the degener- 
acy is lifted, and each multiplet acquires a different en- 
ergy shift denoted by AE S and AE A , respectively (see 
Fig. |6|). We can reproduce this spectrum by consider- 
ing an effective coupling between the spin of the dressed 
spin at site-0, and the spin of the rid electrons on site 
1. If E s lies (above)below E A , the effective coupling is 
(anti)ferromagnetic. 

Thus, if the coupling between the effective spin at the 
impurity site and that of the electrons on site 1 is fer- 
romagnetic we know, from the scaling analysis at weak 
coupling, that the perturbation is irrelevant, and the low 
energy physics is described by the strong coupling fixed 
point. That is, an underscreened, completely symmetric, 
effective impurity weakly coupled to a gas of free elec- 
trons with a phase shift indicating that there are already 
(iV — q) electrons screening the original impurity. In the 
completely antisymmetric case (2S = 1), the phase shift 
corresponds to the unitary limit, S = n/2, for SU (2), and 




Eg 



(a) ferromagnetic effective coupling case 




(b) antiferromagnetic effective coupling case 

FIG. 6: Second order perturbation theory energy shift of the 
strong coupling ground state in the cases of a ferromagnetic 
(Fig. a) and antiferromagnetic (Fig.b) effective coupling 



is a function!^ of q/N for SU(N), reaching the unitary 
limit for q = N/2 (see appendix O). 

If, on the contrary, the effective coupling is antiferro- 
magnetic, the perturbation is relevant, the strong cou- 
pling fixed point is unstable and the low-energy physics 
of the model corresponds to some intermediate coupling 
fixed point, to be identified. The authors in Refs. |24| and 
p5| , have claimed that the flow away from the strong cou- 
pling fixed point is characterized by a two-stage quench- 
ing of the impurity. After the impurity has been partially 
screened by the (N — q) conduction electrons on site 0, 
(N — 1) additional conduction electrons tend to accumu- 
late around the impurity, leading to further screening. 

In this section we explicitely calculate the effects of 
hopping on the strong coupling fixed point to the low- 
est order in perturbation theory, that is, second order 
in t. We will consider the case with an arbitrary num- 
ber nd of conduction electrons in site 1 generalizing the 
case nd = 1 considered in Ref. 25 This will allow us 



to understand the origin of the instability of the strong 
coupling fixed point and eventually enlighten the ques- 
tion about the nature of the Kondo screening : either 
one-stage or two-stage Kondo effect depending on which 
regime is considered. 

Before turning on the hopping term, let us consider 
ground states of the form \GS,nd) — \GS)o\ n d)i, with 
nd electrons on site 1. According to the SU(N) symme- 
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of \GS : rid) s and \GS,rid) A separately, without mixing 
the states. We will thus denote the shifts by AEq and 
AE A , respectively. 

The hopping term is of the form 



Hh = H i + H 2 



t^Tcldn 



+tVdi 



(Hj 



H 2 



FIG. 7: Strong coupling ground state in the presence of n d 
conduction electrons on site 1 coupled either symmetrically 
or antisymmetrically to the dressed impurity on site 



try there are two possible configurations, depending on 
whether the n d electrons are coupled symmetrically or 
antisymmetrically to the composite on site 0. This cor- 
responds to the Clebsch-Gordan series [2S — 1] <£> [rid] — > 
[2S, l 7ld ^ 1 ] e [2S - l,l n *]. We denote the states by 
\GS, n d ) s , and \GS, n d ) A , respectively (Fig. 0). 

The SU(N) symmetry is preserved by the hopping. 
That means that the perturbation will shift the energies 



where S a creates an electron on site 1. We can distinguish 
two types of processes, corresponding to different inter- 
mediate states. The first type, which we denote process 
1, corresponds to an electron hopping from site 1 into 
site first, probing excited states \GS + 1} S ' A , and then 
hopping back to site 1. The indices S, A correspond to 
the two possible intermediate states depending whether 
the conduction electron which hops to site is symmet- 
rically or antisymmetrically correlated with the dressed 
impurity as we will see in details in next section. The 
contribution of the process 1 to the energy shift is the 
following 



2 ^y, (GS, n d \dlcp\GS +l,n d ~ 1)"(GS + l,n d - l\cjd a \GS, n d ) 
1 (E - Ei) 

with i = S, A (see Appendix [^) . 

In process 2, the electron hops from site to site 1 first and then back to site 0, probing \GS — 1), leading to a 
contribution to the energy shift of the form 

2 (GS, n d \cldp\GS -l,n d + 1) (GS - 1, n d + l\djc a \GS, n d ) 



r 



Hence, the energy shifts for the symmetric and antisym- 
metric configurations are, respectively, 



En 



El 



+ 



Mf 



En 



E A 



Mi 



En — El 



(J) 



determines the sign of the effective interaction and the 
stability of the strong coupling fixed point. 



A. Process 1, symmetric configuration 



AEn 



M A 



M 2 A 



En 



E A 



Eq — Ei 



(8) 



where the expressions in the denominators, (Eq — Ef) — 
-AEf, (E a - Ef) = -AE A , and (Eq - E 2 ) = -AE 2 
measuring the energy of the excited states compared to 
the energy of the ground state are given in Table P The 
matrix elements, M, will be introduced below as we will 
study the contribution of each process. The energy dif- 
ference between the two states, 



AEl 



Mf 



A' 



Mf - Mi 4 



A W A — 

~ EQ-Ef ' Ek-Ef 



M 2 S 



M A 



Eq — E 2 



(9) 



We consider first the case where the n d electrons in the 
site-1 are coupled to the site-0 state in the most sym- 
metric configuration. In the shorthand notation that we 
use for the Young Tableaux, it corresponds to the state 
[2S - 1] <g> -► [25, l""*" 1 ]. Here, as opposed to the 
case ri d — lE3, the hamiltonian transforms the ground 
state 



\GS,n d ) s = (d{dl--dl)\GS) , 

into a linear combination of two excited states 
l) s , with energy E$, and \GS 



(10) 

\GS + 
s , with 



l,nd — 1) , with energy E$, and \GS + l,n d — 1, 
energy Ea depending on whether the additional conduc 
tion electron in site is coupled symmetrically or an- 
tisymmetrically to the dressed impurity. The state ob 
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\GS,n d )- 



|GS+l,n<j-l>- 



\GS+l,n d -l)- 



cring hopping of a single electron. From here, we obtain 
the following matrix elements: 



| 5 (G5 + l 1 n < ,-l|ffi|(3'5,n <l ; 



S\2 



{2S + n d -l\ ( 2S + q- 1 



25 



25 + TV - 1 



(12) 



FIG. 8: When rid conduction electrons are coupled symmet- 
rically to the dressed impurity at the origin, the term c^do-, 
acting on the ground state, generates a linear combination of 
two excited states with an additional conduction electron at 
the origin. 



Mf = \ s {GS + \,n d -\\R x \GS,n d ) s \* 
'25- 1A f q-l\ 



t\n d -l) 



25 



N- 1 



(13) 



tained by acting with S a d a on the ground state defined 
by Eq. ^ has to be computed explicitely, and the result 
written as a linear combination of the excited states, Fig. 
|§|. The latter are obtained by coupling the site-0 states 
with (n d — 1) conduction electrons on site 1. The explicit 
expressions for these states are given in Appendix and 
they lead to the result 



£4<k) \GS, 



n d ) 



= 'Ml 2S+ ^~ 1 \GS + l,n d -l) s 



+ Ay/n d - 1 



25-1 
25 



\GS + l,n d -l) b . (11) 



Where the normalization coefficients^!, £1 = \j 2S+N-h 
and A = J -^Py are independent of n d , as we are consid- 



We see right away that M 1 is proportional to (n d — 1), 
and vanishes for n d = 1, whereas M-f depends noticeably 
on n d only for 25 <C n d < N. 



B. Process 1, antisymmetric configuration 

Next, we consider the case where the electrons on site 
1 are coupled to the effective spin on site according to 
[25 - 1] ® [l" d ] -> [25 - 1, l nd ]. In the previous section 
it was easy to write down the strong coupling ground 
state by just putting together the effective spin and the 
n d electrons in the highest weight state possible, to ob- 
tain Eq. (|l0|). Here we have to work out the necessary 
Clebsch-Gordan coef ficien ts. We present some of these 
coefficients in Tabic VII of Appendix The explicit 
form of the ground state is 



\GS, n d )^ bc v 



V/25" 



, / n d + l n d + l n d + l \ 

=_= V2S^ T( J] dl)\GS) aa + II 4)I G5 W > 



(14) 



in the notation y\ — a. As in the n d = 1 case, the hop- 
ping term transforms the state defined by Eq. |lj into 
a state proportional to a given strong-coupling excited 
state (Fig. g), In order to obtain the corresponding ma- 
trix element, we have computed explicitely 



T, C l**J \GS, 



n d ) 



cx \GS + l,n d -l) 



and then we have normalized the resulting state. The 
details can be found in Appendix We have 



E c -M \GS,n d ) 



A yjn~ d \GS + l,n d — 1) 



and the matrix element 

Mf = \ A (GS + l,n d -l\H x \GS,n d ) 

.2 ( 1 - 1 

= t n d 



A\1 



N - 1 



(15) 



Notice the dependence on n d , and the fact that the ma- 
trix element does not depend on 25. Combining together 

M 1 and Mj 4 as it appears in Eq. (||), we have 

Ml - Mi = f 23 ^- 1 ) • (16) 



25 



N - 1 



This is a term with the same n d dependence as Mf but 
with the opposite sign. 
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FIG. 9: When n<j conduction electrons are coupled antisym- 
metrically to the dressed impurity at the origin, the term 
cjyda, acting on the ground state, generates state proportional 
to a given excited state, with additional conduction electron 
at the origin. 
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FIG. 10: The term die, acting on the ground state with 
rid electrons on site 1, coupled symmetrically, (a), or anti- 
symmetrically, (b), to the impurity. The result is a state 
proportional to an excited state with one less electron at the 
impurity site, coupled to (n<j+l) electrons symmetrically, (a), 
or antisymmetrically (b). 



C. Process 2, the trick 

Both in the symmetric and antisymmetric configura- 
tions in process 2 there is a one-to-one correspondence 
between the state obtained by the action of H2 , and the 
excited state with the same symmetrycJ (Fig. [to|). We 
write, for completeness, the Young Tableaux associated 
with these processes The evaluation of the two remaining 
matrix elements, M% and M A , associated with process 2 
is simplified by using the following trick connecting the 
matrix elements of processes 1 and 2, respectively. On 
the one hand, we have in the process 1 

Mf + Mf = t 2 Y, S {GS,n d \dl, ClJ ,cld a \GS,n d ) s 

(7 (T 1 



where SM S = i 2 Xw S (c^d^d^Cc') 3 and the expec- 
tation value 3 () 3 should be considered over the ground 
state, \GS,n d ) ■ On the other hand, the following prop- 
erty holds for the matrix elements of process 2 

Mi = t 2 Y, S (GS,n d \cl,d a/ dU.\GS,n d ) s 

= t 2 J2 s ( c l(^-dUM s 

= t 2 J2 S (c±c„) s -5M s , 

a 

where in the last line, we have exchanged the dummy 
variables a and a' . The same relation exists for the anti- 
symmetric configuration (with expectation values taken 
on \GS,n d ) A ) 

Mi 4 = t 2 J2 A (€da) A -6M A , 
M A = t 2 J2\U«) A -SM A . 

a 

Thus, using the trick, we obtain 

Mi = t 2 (n c -n d ) + (Mi + M?) 



(18) 



N - 1 J \2S + N - 1 



M A = t 2 (n c -n d ) + M( 
= t 2 {N-(n d + l)) 



N-q 
N - 1 



Mo - M. 



A 



Mf + Mf - Mi 4 
2 {2S + n d - 1\ fX 



\2S + N - 1 



N - 1 



(19) 



IV. DISCUSSION 



Once we have computed all the matrix elements, we 
can evaluate the difference in energy shifts between 
the symmetric and the antisymmetric configurations, 
(AEq — AE A ), induced by the perturbation Hh, accord- 
ing to Equation <^). We can then compare the result 
to an effective spin-spin interaction, with coupling J e //, 
between the spin So of the dressed impurity at site 0, 
and the spin Si of the n d electrons at site 1. If the sym- 
metric configuration is the lowest in energy (Fig. [|(a)), 
it means that the effective coupling is ferromagnetic and 
the perturbation is irrelevant. Thus, the strong coupling 
fixed point is stable and allows to describe the low-energy 
behavior of the system. If, on the contrary, it is the anti- 
symmetric configuration the lowest in energy (Fig. ^(b)), 
the effective coupling is antiferromagnetic, and we know 
from scaling arguments that the perturbation is relevant. 
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The strong coupling fixed point is then unstable with 
respect to hopping, and the behavior of the system is 
described by an intermediate coupling fixed point. In- 



corporating the expressions of the matrix elements, Eos. 
( |l4^6| , pj| ) , and those of the excitation energies (Table y) 
into Eq. @, one finds 



fit 2 

AE* - AEq = -(25 + n d -l) I — 

2S + q-l 



+ 



25(25 + N - 1)(25 + N - q - (25 + q - 1)/N) 
(N-q) 



(N — 1)(25 + N — l)(q +(2S + q- 1)/N) 2S(N - l)(N -q-(2S + q- 1)/N) 



r 



(20) 



As the bosonic part increases, the energy difference be- 
comes smaller. The dependence of the r.h.s. of Eq. ( p(j ) 
with (2S + nd— 1) is linear and appears factored out. The 
effect of rid on (AEq — AEq) is weak as long as rid <C 25. 




FIG. 11: Effective coupling, J e ff, as a function of q/N, for 
different values of 2S, in the large-N limit. 

We now calculate explicitely the interaction energy of 
two spins, in the representations [25 — 1] and [l™ d ], re- 
spectively, and identify the effective coupling J e ff from 

^//E( s o 2s -V(sr ] ) A , 

A 

where A = l,N 2 — 1. The spectrum consists of two 
multiplets, according to the Clebsh-Gordan Series [25 — 
l](g)[l" d ]^ [25, l ,ld_1 ] s © [25— 1, Y Ld ] A . We have added 
superindices to indicate which is the symmetric state and 
which is the antisymmetric one. As before, the energy is 



given in terms of Casimir operators, Eq. 
energy difference between the states is 



and the 



AE [2S.V 



l,l n d] 



J, 



e// [C 2 ([25,l" d - 1 ])-C 2 ([25- 



2 

J eff 

4 

J eff 



1,1 



(25 + n d -l) [y'ps.i-a-i] 



Y' 



[25-1,1* 



(2S + n d -l) 



where we have used the results of Table ffl. Since both 
states have Young Tableaux with the same number of 
boxes, Q e ff = 25 + 7id— 1, the energy difference depends 
only on the second constraint (B7), 34// = QeffY' (see 
Appendix |b|) . As a consequence, the dependence on rid 
is factorized in the same way as in Eq. (20). Thus, we 



can identify the effective coupling, J e // with the term in 
curly brackets in Eq. (20) times 4t 2 /J. In Fig. [ll] we 
plot this energy difference for several values of 5. 

Until now, we have presented results for arbitrary N. 
In the following, we will focus on the limit N — * oo with 
2S/N and q/N finite. In this limit the expressions ob- 
tained are simplified, and it is easier to present the main 
features of the model. Thus, in the large-N limit, we have 

•4*2 

. In. 



J eff ~ ~ 



f N(N — 2q) 



1 



V q(N-q) J (25 + N- q) 



.(21) 



J e ff is of order (1/7V 2 ). The result expressed in Eq. ( |2l| ) 
means that the value of the effective coupling J e // is in- 
dependent of the number of conduction electrons, rid, on 
site 1 and coincides with the result obtained in Ref. |25| 
in the case n d — 1 . This is due to the cancellation of the 
(25 + rid — 1) factor that we have mentioned before. The 
coupling remains ferromagnetic as long as q < N/2, as 
can be seen in Fig. 12 and by inspection of the numer- 
ator in the r.h.s. of Eq. (pl|). As soon as q > N/2, the 
strong coupling fixed point becomes unstable. 

The case q = N/2 requires particular attention, since 
the leading contribution to J e tf vanishes. Taking into ac- 
count the whole expression for the effective coupling, we 
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UNSTABLE STRONG COUPLING 
FIXED POINT 

ATTRACTIVE 



STABLE STRONG COUPLING 
FIXED POINT 

REPULSIVE 




0.1 



2S/N 



FIG. 13: Phase diagram of the model (for 2S/N finite), as 
a function of the impurity parameters, 25*, and q. As soon 
as q > N/2 the strong coupling fixed point becomes unsta- 
ble. For q — N/2, the strong coupling fixed point remains 
stable only for moderate values of 2S/N (short line ending in 
a point). 



FIG. 12: Proportionality factor A = N(N — 2q) /q(N — q) as a 
function of q/N . A determines both the sign of J e // and the 
dependence of AE A on n^. 



find that the strong coupling fixed point for an impurity 
with q = N/2 is stable as long as the bosonic parameter, 
S is smaller than a critical value 



S* 




2N 



N- 1 

In the large-N limit we have 

'V2-l\ 4 + \/2 



(N 



S* 



N ■ 



0(l/N) 



N 

To 



We see that the strong coupling fixed point at q = N/2 
becomes unstable already for moderate values of S (Fig. 

It is worth studying the energy shift for each state 
separately. It is easy to show that 



ae£ 



2_e_ 
T 



1 



N — 1 J \ N-q-Q/N 
//,, \ ( N-q_\ 



1 - 



N - 1 



AEg = AE$ + (25 + n d - 1) 



Q/N J 



J, 



eff 



,(22) 



(23) 



Since J e // is 0(1/N 2 ) in the large-N limit, the energy 
difference (AEfi - AE$) is 0(1/N), and both energy 
levels (AE§ and AE^) have the same leading term in 
0(1). Notice that when we consider the large-N limit 
in which Q/N is finite, the antisymmetric energy shift, 
AEq is almost independent of 2S. The dependence of 



AEq on S is mainly contained in the J e ff term, which 
will eventually determine the stability of the system. 

The most important property of Eq. ( f22[) is its behav- 
ior as a function of rid,. To leading order in 1/N, we can 
write 



N-g \ _ (ru\ ( N(N-2q) 
q ) \n)\ q(N-q) 



(24) 



Notice that the factor multiplying n^/N in the r.h.s. of 
Eq. ( |24| ) also appears in J e ff (see Eq. (^l|)). This result 
has the immediate following physical consequence. The 
change of sign of J e // -and hence of the stability of the 
strong coupling fixed point- is directly connected to the 
change in the behavior of AEq ~ AEq with nd- In 
particular, when J e jf = 



9t 2 

AE${q = N/2) = -— , 

for any value of nd- The value of n d that minimizes the 
energy given by Eq. ( |24| ) at arbitrary q will depend on 
the sign of (N — 2q). In the regime where the strong 
coupling fixed point is stable, q/N < 1/2, J e ff < 0, 
the lowest energy corresponds to — 1, whereas for 
q/N > 1/2, J e ff > 0, the energy expressed in Eq. 
is minimized for nd = (N — 1). We have plotted AEq 
in Fig. 14 (compare to Fig. |ll|). The shaded region 
corresponds to the possible values of AEq for the whole 
range of nd, bounded by the limiting cases, nd — 1, and 
nd = (N — 1). In the region where the strong coupling 
fixed point is stable, the electrons on site 1 are repelled 
by the effective impurity (repulsive regime), whereas in 
the unstable, q/N > 1/2 regime, the lowest energy in 
second order perturbation theory corresponds to a state 
where the effective impurity attracts as many electrons 
as possible on site 1 (attractive regime). This is pre- 
cisely the mechanism behind the two-stage quenching. 
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< 




FIG. 14: Leading order term in the energy shift, AEq ~ 
AEg , as a function of g/iV, for 1 < rid < (N — 1) (shaded 
region), and in the limiting cases n d = 1 (dashed line), and 
rij = (iV — 1) (straight line). Notice that the value at q/N = 
1/2 is equal to —2t 2 /J, for any n d . Note that for q/N < 0.5, 
the energy is minimized for n d = 1 while for q/N > 0.5, the 
minimization is obtained for rid = (N — 1) 



The accumulation of electrons on site 1 is not related 
to Jeff which is independent of rid, but results from the 
dependence of AEq ~ AEq with rid. Moreover, as we 
show in Appendix ^], AEq coincides with the energy 
shift for a fermionic impurity (completely antisymmetric 
representation), in the large- N limit. When the impu- 
rity is fermionic, there is no degeneracy of the strong 
coupling fixed point, which is always stable, and the low- 
est order perturbation theory just shifts the ground state 
energy. Nevertheless, there are two regimes, repulsive 
and attractive, depending on the value of q, and char- 
acterized by the value of rid that minimizes the energy. 
This behavior is a consequence of the particle-hole sym- 
metry in the fermionic case, given by the transformations 
q — ► (N - q) and n d -> (N - n d ) (cf. Appendix Oj. The 
behavior of a fermionic impurity with q is the same as 
in the case {N — q), if we reinterpret the electrons as 
holes and the impurity as made out of holes. Therefore, 
if rid = 1 minimizes the energy for q < N/2 (electron 
repulsion), then the energy for a hole impurity, made out 
of (N — q) fermions, is minimized by the state that re- 
pels the holes, (N — rid) — 1, implying an attraction of 
electrons, rid = (N — 1). 

The addition of a bosonic component to the impu- 
rity, leading to the formation of a row in the L-shaped 
Young tableau representing the impurity, breaks this 
particle-hole symmetry. Whereas the two regimes de- 
scribed above are still present, due to the fermionic com- 
ponent, the degeneracy of the states due to the bosonic 



-7V+2S-2- Thus, we should expect a resid- 



lnC* 



N-l 
N+2S-2 



and a Curie susceptdbil- 



C^j + \ s _ 2 /T ', with logarithmic corrections' 



component leads to the instability of the strong coupling 
fixed point at the same point as where the the dressed im- 
purity starts attracting the conduction electrons on site 
1. 

We finish by making some remarks on the physi- 
cal properties of the model in the different regimes. 
As is common to all models with an antiferromagnetic 
Kondo coupling, there will be a crossover from weak cou- 
pling above a given Kondo scale, Tk, to a low-energy 
regime. When the strong coupling fixed point is sta- 
ble, we should expect for T -C Tk a weak coupling of 
the effective impurity at site-0 with the rest of the elec- 
trons. The physical properties at low temperature are 
controlled by the degeneracy of the effective impurity, 
d([2S-l]) = Ci' 

ual entropy S l 

ity, x l 

This is the result that we would expect for a purely sym- 
metric impurity. The difference with respect to the case 
at hand is that in the L-shaped impurity model only 
(N — q) electrons are allowed at the origin, instead of 
(N — 1). Thus, we would expect to find different results 
for quantities that involve the scattering phase shift of 
electrons off the effective impurity (see Appendix 

In the q > N/2, we do not have access to the intermedi- 
ate coupling fixed point that determines the low-energy 
behavior. Nevertheless, it is reasonable to think that 
there would be a magnetic contribution to the entropy, 
and a Curie-like contribution to the susceptibility, since 
the impurity remains unscreened. This behavior is differ- 
ent from that of the multichannel Kondo model, which 
is also characterized by an intermediate coupling fixed 
point, but where the impurit^-jpagnetic degrees of free- 
dom are completely quenched^. The degeneracy of the 
true ground state is an open question, but we can assume 
that the entropy wilLbe smaller than that of the strongly 
coupled fixed pointo. It is in the scattering properties 
that we might be able to see the anomalous features of 
this new fixed point more clearly. 



V. CONCLUSIONS 

In this paper, we have studied the SU(N), single- 
channel Kondo model, with a general impurity spin, in- 
volving both bosonic and fermionic degrees of freedom 
(corresponding, respectively, to the horizontal and verti- 
cal directions in a L-shaped Young tableau) . This model 
shows a transition from a strong coupled fixed point to 
an intermediate coupling fixed point, when the amount 
of fermionic degrees of freedom, q, becomes larger than 
N/2. We have identified the origin of this instability as 
due to the change, from repulsive to attractive, of the 
effective interaction between the dressed impurity and 
the conduction electrons in the neighboring sites. This 
change is already present in the purely fermionic case, 
where it happens at the particle- hole symmetry point, 
q = N/2. 
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These results shed light on the nature of the Kondo 
screening. Whereas the screening is restricted to a one- 
stage process in the strong-coupling regime, involving 
mainly the fermionic part of the impurity, the accumula- 
tion of conduction electrons on neighboring sites in the 
attractive regime gives rise to a two-stage Kondo effect. 
The latter case corresponds to a regime where the low- 
energy physics is controlled by an intermediate-coupling 
fixed point. Obviously, the interesting open problem now 
is to understand the physics associated with the this 
new intermediate coupling fixed point. This issue might 
have important future applications for the lattice prob- 
lem, with potential consequences for the understanding 
of non-Fermi liquid behaviour observed in heavy-fermion 
systems. 
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APPENDIX A: COMPOSITION OF THREE 
FUNDAMENTAL REPRESENTATIONS OF SU(3) 

Before dealing with the general problem of construct- 
ing the high est weight impurity states in SU (N) in Ap- 
pendix B 1 , we will write down in detail all the three- 
particle states with SU(3) symmetryEa. This will al- 
low us to see how the states constructed with different 
numbers of bosons and fermions can be basis for repre- 
sentations with the same Young tableau. We will also 
see the role of the SU (3\3) and 5£7(1|1) supersymmetry 
groups induced by the realization in terms of bosons and 
fermions. 

The direct product 3 <£> 3 <£> 3 of three fundamental rep- 
resentations of SU (3) gives the following Clebsh-Gordan 
series 



3® 3® 3 = (3® 6) ©(3® 3) = lOeS^S 2 ©!. , (Al) 

where we identify each representation by its dimension. 
In terms of Young tableaux, we have 



The representations 6 and 3 result from the composition 
of two fundamental representations 



3(g3 = 6©3 . 



(A2) 



In addition to Young Tableaux, we can use weight dia- 
grams to describe the states in the representation, Fig. 
53 . In SU(3), we associate a triangle to the fundamental 
representation, 3. Each vertex corresponds to a partic- 
ular state of the multiplet, and the different states are 
related by the action of the lowering operators. In Fig. 



IbV - 



T 



la? 




IcV 



FIG. 15: Weight diagram for the fundamental representation 
of SU(3), with the states and the relevant lowering operators. 

|l^we include the weight diagrams associated with Eq. 

The representation 6 is completely symmetric. Its 
states are realized in terms of Schwinger bosons. For 
instance its highest weight state can be written as 



1 



(bl) 2 \0) 



Here, and in the following, the values of SU(3) spin are 
denoted by a, 6, and c. Likewise, the representation 3 is 
completely antisymmetric, and its states are more con- 
veniently expressed in terms of fermions. For the highest 
weight state we have 

|a6)f =/t/t|0) . 

There is another way of realize both 6 and 3 using one bo- 
son and a fermion. Being symmetric, the highest weight 
of 6 is easy to write, since 



is already symmetrized. To obtain the highest weight of 
3, we first have to find a state with the same quantum 



numbers in 6, by acting with T ab = (Jlf a + b\b a ) on 



to get 



/M)|o) 



(bb) 6 (ab)6 (aa) 6 (ab)3 




fee) 6 

: Weight diagram for the Clebsch-Gordan series of the product of two fundamental representations of SU(3) 
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and find a state orthogonal to \ab)J 

\ab)f = -L(ftbi-ftbl)\0). 

This process is described in Fig. |6[ It is easy to see that 
the states with the subindex / are related to those with 
the subindex b by the SU supersymmctric operator 

= bo that <?|(.--))/=V2 |(- ••)>&• 

Consider now the three-particle states. The easiest 

state to write is the highest weight state in the most 

symmetric representation, 10 (cf. Fig. |l7| ) 

\aaa) 10 = \a) 3 \aa) 6 . 

It can be e asil y expressed in terms of bosons in agreement 
with Eq. (pl|) 



)l° = 1 (6t)3| 0> 



Alternatively, we can use a realization with two bosons 
(from 6) and a fermion (from 3) 

,f = -!=/M) 2 |o> . 

V 2! 



Other states of the representation are obtained by the 
repeated action of lowering operators. For instance, 

\aab) w = -^(V2\a) 3 \abf + \bf\aa) 6 ) , (A3) 
v3 



which leads to 



\aab)l° = -L(bifbl\0) 



ib)f=±(2Mbi+ft(btr) 



The octets, 8 are mixed symmetry representations that 
have to be built by a combination of fermions and bosons. 
The Clebsch-Gordan series (A.1) indicates that 8 1 , built 
from the product 3 <S> 6 is naturally realized by states 
with one fermion and two bosons, Fig. [l7|, whereas 8 2 is 
realized by the product of one boson and two fermions 
(3 (g) 3), Fig. |l|. The highest weight state, \aab) 8 \ of 



the octet 8 1 is orthogonal to \aab) 10 defined in Eq. (A3) 



\aab)f = ±(\a) 3 \ab) 6 - V2\b) 3 \aa} 6 ) 

= ±bUtibt-fibl)\0) , 

in agree ment with the general expression of the state ipb 
given in (B3). 

As usual, the other states of the octet are built by 
the repeated action of generators = (flf a + b\b a ), 



U- = (fife + b\b c ), and V- = (/J/ c + of 5(7(3). 

For instance 



) 8± 
lb 



1 



\8 X 



V2 J 



= -j=(\a) s \bc) 6 + \b) 3 \ac) 6 -2\c) 3 \ab) 6 ) 

= ^(flblbl + ftbibt-2Mbt)\Q) ■ 

This state is degenerate, since there is another state in 
the multiplet with the same quantum numbers. In or- 
der to find this last state, \acb) 8 , we have to combine 
the action of lowering operators with orthogonality with 
respect to \abc) 8 . Acting with U~ on |a6£>) 8 leads to 
a state that is not orthogonal to |a&c) 8 . Therefore, we 
write 

U-\abb) &1 = V2(a \abc} 81 +(3 \acb) sl ) . 

There are two ways of reaching states with quantum num- 
bers {abc} starting from \aab) s . We use this and the fact 
that [T+, U~] = 0, (T+ = (T")t), to derive a 

2 a = V2 8l (abc\U~\abb) al = 8l (aac\T + U-T-\aab) 81 
= 81 (aab\U + T + U-T-\aab) &1 



8 1 



(aab\T + T~U + U-\aab) 



8 1 



1 



Hence one deduces a = 1/2 and = V3/2 and one gets 
the expression of the last state \acb) 8 of the octet 



\acb)f = -^(|a) 3 |6 C ) 6 -|6) 3 |a C ) 6 ) 
= ^(/+6i6t-/t 6 t 6 t)| 0> . 

To summarize, the rule is to describe the highest 
weight state of the octet \aab)b using the bosonic rep- 
resentation, ij)b, and hence to derive the other states of 
the octet by the above construction allowing to recover 
the Clebsch-Gordan coefficients involved in the spin com- 
position related to the direct product 3 <g> 6. 

The highest weight state of the octet 8 2 is also the 
highest weight state of 3 (g> 3 

\aab)f = \a) 3 \abf = bUlft\Q) , 

i n a greement with the general expression of the state ifif, 
(B4). The construction of the other states follow the 
same lines as in the case of the 8 1 octet. For instance, 



1 



\abc)f 



V2 



(\af\bcf-\bf\ca) 3 ) 



^{bVtfi-bififtm , 



\acb)f 



1 



(\a) 3 \bcf + \bf\caf-2\c) 3 \abf) 
(blfif±+btf±ft-2btftft)\0) . 
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Once again, the two basis of states corresponding to the 
same Young tableau, are related by the SU(1\1) opera- 
tors, 9\ and 9. 

Finally, the singlet state, \cibc) 1 , is built by orthogo- 
nality with the states | abc) 8 ' , and \acb) s , from the octet 
in the product 3 <g> 3 (Fig. Q , 



1 



\cxbc) 1 = -j=(\a) 3 \b 



\bf\ca) 



lab)' 



(A4) 



in agreement with the expression of the states in the com- 
pletely antisymmetric representation of the spin ( |B2| ). 
The simplest way to realize this state is with three 
fermions. Then 



\abc}) = ftfl ft |0> 



(A5) 



But it can also be written with one boson and two 



fermions, either by acting with 9 on (A5) or by substi- 
tuting on (A.4) 



\abc)l = -^{btftft + blftfl + blftft) . 

We have shown how to construct the states for dif- 
ferent representations of SU(3). These results can be 
summarized as tables of Clebsch-Gordan coefficients like 



the ones that we presented in Tables IV -VI] for SU(N). 

Let us make some comments about the number of 
states. The direct product of three fundamental repre- 
sentations of SU(3) generates a space of dimension 27, 
which breaks down as a direct sum of irreducible repre- 
sentations according to the Clebsch-Gordan series (Al). 
By considering all the realizations of these states in terms 
of bosons and fermions subject only to the constraint 
Q = rif + rib — 3, we are working on the higher dimen- 
sional space of a representation of SU (3|3) with a total of 
38 states, as represented in Fig [l9] .The figure also reports 
the relation between these states with the introduction 
of an additional supersymmetric operator, Z~ r = ft b a 
, acting on the highest weight states. The states 10 and 
10' respectively are identical (as well as 1 and 1') as far 
as the SU(3) symmetry is concerned. This is not the 
case of the two states 8 1 and 8 2 which correspond to dif- 
ferent spin representations even if the associated Young 
tableaux are the same. Altogether, one recovers the ex- 
pected total of 27 different states. 



m = 3, n/ = 10 | b \ b \ b \ 

z- 1 



n b = 2, n f = 1 8 1 



LfJ 



10 ' I / 1 b | b | 



0+ 



rij, = 1, nj = 2 1' 



b 

L 
L 



0t 



n b — 0, rif = 3 



in 



JL 
L 



FIG. 19: States with Q — (rif, + rif) grouped according to 
SU(3) representations 



1. Impurity State 

Before studying the general case of a L-shaped Young 
tableau representation of spin, we will consider the two 
limiting cases of a completely symmetric (bosonic) and 
antisymmetric (fcrmionic) representations of the spin in 
SU(N). 

The case of a completely symmetric representation of 
the spin is equivalent to a system of 25" identical particles 
symmetric under the permutation of two of them. The 
associated Young tableau is made of a single line of 25* 
boxes 



2.S' 



Pes 2S 



expressed in shorthand notation as [2S] . Associated with 
the Young tableau, there is a symmetrizer operator made 
out of the sum of all the permutations of 2S elements. 
It is convenient to use an explicit representation of the 
localized spin in terms of N species of Schwinger bosons 
b a , (a — a, 6, ..^r^) subject to the constraint 



APPENDIX B: IMPURITY STATES AND 
LOW-LYING STATES IN THE 
STRONG-COUPLING LIMIT 

We present in this appendix the explicit form of the 
highest weight (spin) states for the impurity and for the 
low- lying states in the J — > oo limit. We discuss in detail 
the use of Young Tableaux both to describe the sym- 
metry properties of the states and to study de SU(N) 
generalization of the composition of several spins. 



a 

The (N 2 — 1) components of the spin operator can be 
represented as S A — Y^ a d ^a T apbf3 1 while the highest 
weight state (the analog of the state with the largest 
value of S z in 5/7(2)), can be written as 

\{af s )^ = ^={blf s \Q) , (Bl) 
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where |0) denotes the vacuum state for the bosons. 
Other states of the representation can be obtained from 
this one by the repeated action of lowering operators, 
taking advantage of the underlying SU (2) subalgebras 
within SU(N). Take, for instance, the (25 + 1) states 
{\(a) x (b) v )^ 2S ^} with x + y = 25. They transform as 
a regular, SU(2), spin-5 multiplet under the action of 
the SU(N) operators T~ b = b\b a , T+ = b\b b , and 
\bb)/2. In particular, 



a b 



(blba 



T- \a 2S )^ =V25 Ka) 23 -^) 



2S] 



or, in terms of bosons, 



V(25-l)! 



i) (95 - 1} 6j|0> 



Note that each index in the set of quantum numbers, 
{a, /?..., p2s}, describing the states of the representation, 
can take N values independently of the rest of the set, and 
that for each set of values there is only one state. The di- 
mension of the representation is thus given by Cjf, 2 s— v 
corresponding to the number of ways of choosing 25 ele- 
ments out of a group of (N + 25 — 1). 

The other limiting case corresponds to a completely 
antisymmetric representation of the spin. It is equivalent 
to the case of q identical particles antisymmetric under 
the permutation of two of them. The associated Young 
tableau is made out of a single column of q boxes 



with respect to these subalgebras. The dimension of the 
representation is given by Cfj. 

An important property of both kinds of representations 
is that the states are non-degenerate. That is, each set 
of allowed quantum numbers completely determine the 
state. This is not the case for mixed symmetry represen- 
tations, as we are about to see. 

Let us now consider the general L-shaped represen- 
tation of spin, Fig. [I], which interpolates between 
the previous two limits. Its dimension can be eas- 
ily obtained using Robinson's formulae^: the result is 

f_2S \ s-i2S /-Y9-1 



r<2S 

Ly AT+2S- 



{ 2S+ q -i J ^n+2S-i C n-i (Table 0). The L-shaped rep- 
resentation is the result of the direct product of a sym- 
metric and an antisymmetric representation. This can be 
done in two, non-equivalent ways (Fig. |20| ) : either as the 
[l 9-1 ] <g> [25] -> [25, l q - x ] © • • • Clebsch-Gordan series, or 
as [l q ] <g> [25 - 1] -> [25, l^ 1 } © • • ■ . The construction of 
the highest weight states for each of the cases is d etaile d 
at the end of this subsection (cf. Eqs. (BE) and ( B12p ) 
and leads to 



E 



p 



Pes„ 



and expressed in shorthand notation as [l 9 ] with q < N. 
Associated with the Young tableau, there is an antisym- 
metrizer operator made out of the sum of all the per- 
mutations of q elements weighted by a Sp = ±1 factor 
as for antisymmetric identical particles. It is convenient 
to use an explicit representation of the localized spin 
in terms of N species of Abrikosov pseudo-fermions f a 
(a = a, 6, ...,rjv) subject to the constraint 



ft fa 



q 



The generators of SU(N) in this realization are S A = 
J2 a g fa T afifl3i an d the highest weight state of the repre- 



a/3 

sentation can be written as 



|o6...r,)[ 1, ]=/t/t.../t| ) 



(B2) 



involving a set {a,b, ...r q } of q different indices. Other 
states of the representation can be obtained from this 
one by the repeated action of lowering operators such as 
T~ b = flfa, taking advantage of the underlying SU(2) 
subalgebras. Note that the states are SU(2) doublets 



FIG. 20: Two ways of obtaining a L-shaped representation out 
of the direct product of a symmetric and an antisymmetric 
SU(N) representations 



06 



and 



\2S 



b...r, 
1 



\[2S , ,1 5 ~ 1 ] 



(63.) 



2S-1 



y/2S + q-ly/(2S-l)\ 



^ = \(a) 2S b...r q ^ 2S ^ 



■A{b\ 4/1-4)10) , (B3) 



(6 f J 



,2S-1 



V(25-l)! 



{flflft-flm , 



(B4) 



where A(- • • ) is the antisymmetrizer. The impurity spin 
operator has a form which is independent of either 25 or 
q and is given by 5 = J^apibaTapbp + fo- a f,f fi ). The 
first constraint has to do with the conservation of the 
number of particles 



Q = (hf + h b ) = (2S + q- 1) 



(B5) 



where (25 + q — 1) denotes the number of boxes in the 
L-shaped Young tableau. In the limiting cases discussed 
previously, once the value of Q is fixed, the representa- 
tion is completely determined. Here, however, it is nec- 
essary to add a second constraint to identify states with 
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the right symmetry. To that end, it is worth noticing 
that the set of states tpb and ipf form a basis for a repre- 
sentation of the larger, supersymmetric group SU(N\N), 
with generators given as linear combinations of the op- 
erators 6+6,3, ftfp, bifp, fib p. Thus, all the L-shaped 
impurities that interpolate between the symmetric and 
the antisymmetric case are related by the supersymmet- 
ric group. As a matter of fact, the constraints that fix 
the SU(N) representation are obtained from operators in 
SU(N\N) diagonal in spin, such as fif, n b , 9 — VJ a 6^/q, 
and 9^ . Consider for instance, the action of 9' on ipb, 

0t (6t)^-i^t / t / t... / tj| 0) 
= (25 + g -l)(6t) 2S - 1 (/t/t/t.../tj |o) . 

where the right hand side corresponds to ipf. This leads 
to the relations 



(B5) is fulfilled, the invariance of the Casimir C2 is en- 
sured provided that the operator y is invariant too. This 
leads to the second constraint 



y = Q(q-2S) 



(B7) 



It is easy to check that the operators 9 and 9* commute 
with Q and y , which implies that the constraints are also 
compatible with the SU supersymmetry. Note that 
this was not the case for the operator y' — fif — hb + 
^[0, 0*<} defined in Eq. 10 of ref. ||. 

The constraints completely determine the representa- 
tion, but they cannot distinguish between the states tpt 
and ipb- The physical properties of the system depend 
only on the Young tableau associated with the Kondo 
impurity, and not on the particular way the representa- 
tion basis is constructed. 



0^ b = ^/2S + q-liP f , 
Ipf = ^2S + q-lfo . 

Notice also that 9'tpf — 9ijjb = 0. The operators 9 and 
9^ relate states that transform under a representation 
of SU(N) given by the same Young Tableau. Together 
with Q, they form the 5/7(1 11) supersymmetric algebraEJ 
{9,9^} = Q. Furthermore, the operators Pb = -^99^ 
and Pf = ^9^9 arc the projectors out of the bosonic 
ipb and the fermionic ipf states, respectively. The states 
■0 / and ipb are the exact analog of the familiar example 
of the formation of the two octets, 8 1 and 8 2 , out of 
the composition of three fundamental representations in 
SU(3). In appendix [A| we construct the states explicitely 
in this example and derive the corresponding Clebsch- 
Gordan coefficients. 

The second constraint is then given by y, a bilinear 
combination of the operators {n/, n;,, since it is a 

consequence of the invariance of the Casimir operator, 
C 2 , (the SU(N) generalization of S 2 = S{S + 1)), which 
for a L-shaped representation is given by 



Q(N-^)-y 



(B6) 



Here, y — Q{hf — fib) + [9, 0']. Once the first constraint, 



We finish this section by describing in detail the con- 
struction of the relevant states of the impurity multiplet. 

The direct product of irreducible representations of 
SU(N) decomposes into a direct sum of irreducible rep- 
resentations (Clebsch-Gordan series). A well-known ex- 
ample is the addition of angular momentum in SU(2). 
In order to find the states in the new basis (leading to 
the Clebsch-Gordan coefficients (CGC)), we follow a sim- 
ilar procedure to the one used for angular momentum 
(see also Appendix |a|). That is, we first identify the 
highest weight state of the product of representations, 
which is always non-degenerate, with the highest weight 
state of the most symmetric of the representations in the 
Clebsch-Gordan series. Then, we use lowering operators 
(particular combinations of the generators of SU(N)) to 
generate other states in the same representation. For ar- 
bitrary N there might be more than one state with the 
same quantum numbers in the same representation (as 
in the SU(3) octet, 8 1 and 8 2 ), and we will have to find 
orthogonal combinations. Finally, we find states in the 
next representation by looking for additional orthogonal 
states with the same quantum numbers, and acting on 
them with lowering operators. 

The direct product of a symmetric and an antisymmet- 
ric representation of SU(N) can be expressed as a direct 
sum of two L-shaped representations, according to the 
Clebsh-Gordan series 



[25-1] 

[2S-1] 



[25, V'- 1 } 



[25- 



(9 + 1) 
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with dimensions 



/-as -i r<q _ 25 

^N+2S-2 2S + g _ 1 



r<2S 

L/ N+2S- 



n q-\ 
1 °JV-1 



25-1 



25 



C 



2S— 1 

7V+25-2 ^N-l 



We proceed now to write the states in each L-shaped 
representation. The most symmetric (highest weight 
state) of the product [25 - 1] ® [l 9 ] is 

\(ar-if s -V\abc...r q rK 

where 25 particles have the same quantum number a. 
This state is also the highest weight state of the most 
symmetric of the representations, [25, 1 9_1 ], 



It is easy to see that it must be of the form 

Ka) 2S - 1 6cd...7 V • 9+1 ) [2S - 1 < 1,1 
= 1 (V2S~T | (a) 23 - 1 ) t 2S -^\bcd ...r q )^ 

\j2S J rq—l 

(B12) 



\{af s bc...r q )^ iq ^ = \(a) 2S '^ 2S ^\abe...r q )^ . 

(B8) 

Notice that the lowering operators that transform the 
value a into 5 e [b, r q ] affect the [25 — 1] state only, since 
the [l q ] term is completely antisymmetric. Thus, these 
states are also non-degenerate 



x[2S,l q 



\( a y^bc...(sy... lq/ 

= \(a) 2S - 2 5f s -V\abc...r q )^ . (B9) 



2. Ground state 

The explicit form of the highest weight state for the 
ground state corresponds to Eq. (H), 



with 



\A 2S -!) ! 



Other states with 25 of the ay equal to a, can be ob- 
tained from Eq. ( p38| ). They are also nondegenerate. For 
instance the state with the value 7 replaced by r q+ \ is 



a=i2 /3=i,4 



\(a) 2S bed...r q r q+1 (no -f))^^ 

= KafS-^-^abcd. . .r q r q+1 (no 7))^] .(BIO) 

Altogether, there are g states of this type (if we restrict 
the ctj to the range [a, r q+ {\). 

Next, we construct the states in [25, l 9 ^ 1 ] where (25 — 
1) of the ctj are equal to a and the rest of labels are 
different and take values in the range [6, There are 

q linearly independent states of this kind. For instance, 
acting with T~ e on ( B10| ) we get 



T ac\( a 



,2S 



bed . . . r q 

2S-1, 



q r q+1 (no 
2S\(ar^bed...r q r q+1 )^ iq ^ 
V2S~1 \(a) 2S - 2 ^ 2S -%be...r q r q+1 (no j))^ 
\( a) *s-i ) \*s-i] ]bcd ...r 9 r 9¥1 )M . (Bll) 



The orthogonalization of these q states leads to the first 
q rows of Table [\T| of CGC. The last line in the table 
correspond to a state with the same quantum numbers 
((a) 2S ~ 1 bcd . . - r q r q+ i), and orthogonal to all the states 
in [25, l q }. This state is the highest weight of [25-1, 1«]. 



1 =y/(2S + N-l)C q N -- 1 ■ 



Notice the additional term C% \ in the normalization 



factor 7, as compared to Eq. (B3), due to the presence 
of two kinds of fermions f\ and c^7 We adopt the bosonic 
realization of the impurity tpb, which simplifies the calcu- 
lations. We would like to emphasize that all the results 
are independent of the realization chosen, since only the 
form of the Young tableau is relevant to the interaction. 

Other states in the same [25 — 1] multiplet can be ob- 
tained by just acting on the (b\) 2 term. For instance 



T \GS){ a } aa = \/2S - 1 |G5){ a } ab 



(25 - 1) 



1 



bl(b{) 2S - 2 \A) , (B13) 



U |G5){ a } ab = |G5){ a } ac 



1 



bt(b{) 2S ~ 2 \A) , (B14) 
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and 



\GS){a}ab 

\GS){ a } bb 
\GS){a}ac 



- 7 =i= ( 6D=(6t)«-|A) 



1 



\A 2S - 2 ) ! 



L=6j6t (6 t)M-8| A> . (B15) 



V(2S=3)T 



Excited states 



Let us now write the expression of the excited states of 
the strong-coupling fixed point. The states in \GS + l) s 
transform as the completely symmetric representation 
[25]. The highest weight state can be obtained by acting 
with c\ on the ground state, 



\GS + l)f a } oaa — qC^|G?S){ } 00 , 



2S+q-l 
2S+JV-1' 



where the normalization factor, ft — 

appearsil because the additional c-electron has to be an- 
tisymmetrized with respect to the (N — q) electrons al- 
ready present on site-0. Other states in the multiplet can 
be obtained by repea ted action of the lowering operators, 
as in Eqs. (B13,B15). For instance 



\GS + 1 
1 1 



{a}aab 



V2S^Tc{\GS) {a}ab + cl\GS) {a}aa hl6) 



\GS + l)f a } ahc 



1 1 



n V2S 

Cb\ GS ){a.} 



V2S^2cl\GS) {a}bc 
- cl\GS){ a y ab 



States in \GS + l) A transform as [25-1, 1]. The highest 
weight state is \GS + l)f a y aab , and it will be orthogonal 
to the state defined in Eq.(B16). Thus 



\GS + l)f a}aafc 



1 1 



A V2S 



cl\GS) {a}ab -V2S~lc\\GS) {a}o 



wit hcl A = Jjj^. Notice the difference with respect to 
the normalization factor. Only depends on 25. 

Other states in \GS+1) A are obtained in a very similar 
way to the construction of the octets, 8 1 and 8 2 in SU (3) 
(see Appendix . For instance 



IGS+l) 



{a}abc 



-c b \GS) ac 



A y/2S(2S~T 
(2S-l)cl\GS) ab ] , 



V2S^2ci\GS) t 



\GS + 1) 
1 



{a}acb 

l 



A V25 - 1 



c{\GS) bc -V2S^2cl\GS) a[ 



We collect the coefficients for the states with quantum 
numbers {{ajafec} in Table III of Clebsch-Gordan coeffi- 
cients corresponding to the direct product [1] ® [25 — 1] 
in SU (N) , following the notation of ref . ^4|. 

TABLE III: Clebsch-Gordan coefficients for the excited states 
with one additional electron, corresponding to the product 
[1] CS> [2S — 1], and to states with quantum numbers {{a}abc}. 
The normalization factor is l/y/~N, and a is understood 
over each coefficients. The minus sign indicates a negative 
sign in front of the square root, — 





N 


cl\GS) {a}bc 


ct\GS) {a}ac 4\GS) {a}ab 


n\GS + i 


s 

{a}abc 


25 


2S - 2 


1 1 


A\GS + 1 


A 

{a}abc 


2S(2S - 1) 


2S - 2 


1 -(2S-1) 2 


A\GS + 1 


A 


2S - 1 


1 


-(25-2) 



We finish this section with the excited states corre- 
sponding to the multiplet \GS — 1), characterized by one 
less conduction electron than in the ground state. They 
transform as [25— 1, l^ -1 ], and the highest weight state 
is 

1 



\GS-1) = 



(bl) 



2S-1| A / 



with 



l q IjV-l 



i = ^(25 + N - 2)C q N _\ . 

Here, |A') transforms as the JV-dimensional, [l^ -1 ], rep- 
resentation of SU(N). 



APPENDIX C: THE CASE OF THE 
ANTISYMMETRIC IMPURITY 

The ground state for the strong coupling fixed point of 
a Kondo model with a fermionic impurity (completely an- 
tisymmetric representation), \GS)^ 1 is a singlet formed 
with (N — q) conduction electrons at site 0, and has the 
same form as |A), in (|^), but with a different normaliza- 
tion factor 



\gs)^ n ^ ^=A((]j /t)( n . 



(CI) 



P=ia 
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The hopping hamiltonian leads to two types of processes, 
as described earlier, where the intermediate states have 
either one more or one less conduction electron. With 
the help of Table | we can write the excitation energies 

J / N + 1\ , , J /JV + 1\ 

It is clear that AEi and AE2 are related by the particle- 
hole transformation q <-> (JV — g). 

The ground state with n<j electrons on site 1 is 

\GS,n d ) = (dtdl--dt)\GS) . 

Since there is only one intermediate state for each pro- 
cess, we can use the trick described in section [II C, and 
write 



Mi = t 2 



- ^(G5, n d \e\d\,d a e a , \GS, n d ) J , (C2) 



As long as n d > 1, the second term in Eq. (C2) vanishes 
for a 7^ & '. c a acting on \GS,n d ) just counts the number 
of terms where there is a c\. There are C^ r _ 1 such terms. 
Thus, 



J2(GS,n d \cidl,d a c a ,\GS, 



n d ) = n d 



= n d 



( C 9 N-1 

N — q 

JV 



and 



Mi = t 2 n d — 
N 



Finally, following Eq.(p 



M 2 = t 2 (n c - n d ) + Mi = t 2 {N - n d ) 



N-q 
JV 



Mi and Mi are related by the same transformations 
(q <-> JV — q, n d <-> JV — n d ) as the excitation energies. 
That means that the energy shift is invariant under these 
transformations 



AE f = - — 



n d 



N + l \N 



q 



N-n d f N-q 
JV+1 V Q 



In the large-N limit, this result is equivalent to Eq. (24). 
To leading order in 1/JV, the energy shift of the strong 
coupling fixed point is determined by the fermionic com- 
ponent, and the behavior under the particle- hole trans- 
formation. 

We end this appendix by studying the phase shift (5, of 
the conduction electrons scattered off the impurity site. 
This quantity characterizes the impurity contribution to 



the resistivity_p l . At zero temperature and magnetic 
field, we havec3 



p oc sin 



2 6 



(C3) 



The phase shift for antisymmetric impurities in SU (JV) , 
was computed in ref. |3^. In the completely screened case 
it reads 

e 2i5 = _ e -«r(l-2^) _ ( C4 ) 

If we choose the phase shift so that \5\ < f , we have 
tt(£), g<JV/2 



5 = 



-tt( V). 1> N / 2 



(C5) 



The unitary limit, \S\ = 5 is reached in the particle- hole 
symmetric case, q = N/2. We see that this corresponds 
to the point where AEf is independent of n d , indicating 
the change from the attractive to the repulsive regime. 



APPENDIX D: DETAILS OF THE 
CALCULATIONS OF THE MATRIX ELEMENTS 

Here we construct explicitely the excited states that 
are involved in the second order perturbation theory, and 
then we compare them to the action of the hopping term, 
(clda), on the ground state. First, we add a c-electron to 
the ground state, and then we combine it with (n d — 1) 
electrons from site 1. 



1. Symmetric process, \GS, rid) 

The strong coupling excited state \GS + l) 5 is easy 
to compute, since it is the highest weight state in the 
product [2S - 1] <g> [1] -> [25] © [25 - 1, 1]. We have 

|G5 + l)f aa = i C t|G5) ao , 

with the normalization factor fl 2 = 2S+n~-i • Other 
states within the same representation, which transform 
as [25], can be obtained by acting with the corresponding- 
lowering operators. For instance, 

\GS + l) s aab = -±= (V25^T ct|G5) Qb + el\GS) aa ) . 



The highest weight state of the antisymmetric multiplet 
|G5 + 1) A , which transforms as [25 —1,1] is a state or- 
thogonal to |G5 + l)f ab , that is 



\GS+1) 



l j= (c\\GS) ab - V2S~l 4\GS) aa ) 



AV25 
it norm; 

states can be obtained from these three. The results are 



with a different normalization factor, A 2 = jf^j- Other 
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summarized in the Tables [V-VI of Clebsch-Gordan co- 1. The result of the calculations are summarized in the 



efficients. Tables [VI]] . We have 

Next, we have to add the (rid — 1) electrons on site 

I 



\GS + 1, rid - l)oa&— u — ^2g- 



1 / ng n d rid \ 

=_= [V2S (YldDlGS + l^ + J^i-^i J] dl)\GS+l) s aax \ 



1 



25(n4j4|G5) Qa + ^(-ir 1 ( [] 0(V2^T4|G5) a:C3 +4jG5) aQ )] , 

i=2 3=2 »=l,*7 t j 



0^25(25 + ^-1) 
for the symmetric state in the symmetric configuration, and 



|G5+l,n rf -l)f a6 ...„ 



1 



1 / " d . ' id \ 

==[j2(-iy( n 4)^5+1)^ 



A A /25(r^ - 1) 



Et-IR II 4j(4|G5) a ^-V25^I4jG5) aa ) 

j=2 i=l,»^j 



for the antisymmetric state in the symmetric configuration (x\ = a). From here, it is easy to show that the effect of 
{cld a ) on the ground state, is given by Eq. ([Tl]) . 

2. Antisymmetric process, \GS, rid) A 
The action of (cJ.d CT ) on the ground state \GS, rid) A , with rid electrons on site 1 coupled antisymmetrically, produces 



E4<k \GS,n d )$ h ... v = I ! V25^1 Y, (-!)'( II 4)4,|G5) aa 

"d + l [ j 1 rid + l rid + 1 rid + 1 I \ 

+ Ef-^MEf- 1 )'"^ II 4)+ E (-!)'( II 4) \<\Gs) ax A , (di) 



which is proportional to a given strong coupling excited state. Since we can write (Dl) as 

' rid + l rid + 1 

V25 + n d 



' v25+n rf -l 1 Is i=27^ 



+ 1 / j— 1 tid+1 rid + 1 

^(-r 1 E(-!)^( II 4) - £ II 4 4 )]4,ig5)^ 

3=2 \i=2 i=l,ijtj,l 1=3+1 i=l,ijtj,l 



after some algebra, we get 



\n d -l n '' + 1 ! "<i+ 1 



E4^l |G5,n d )l..„ = -!=£=_ E (-!)'{ ( II 4,)(V25^T4 i |G5) QQ -4|G5) aa;! ) 

l-l n d +l 

- £(-l) J ( II 4,)(4jG5) a . 3 -4jG5) c 

j=2 i=l,i^3,l 
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This expression can be written using the antisymmetric states \GS + 1} A , with the help of the following relations 
A(V2S \GS+l)± [Xj - V2S^2 \GS+1)±. XI ) = V25=T(4, \GS) aXj - cl.\GS) axl ) , 



AV2S \GS + l)f = c{\GS) axl - V2S^T 4, \GS) C 



to obtain 



E4 d <x \GS,n d )£ b ... v = 



(_l)nd-iA 



y/(2S + n d -l){2S-l) 



n-d + 1 1 n d + l 

xEH)' +1 ( II 4)^25(25-1) |G5 + l)f a 



i=i 



Finally, since 



EC" 1 )^ II 4j(V25|G5 + l)4. Ki -V2S^2|G5 + l)^ iX ; 

j=2 *=M#j,Z 



X]4^j \GS,n d )i b ... v oc |G5 + l,n d - 1)^...,, 



(D2) 



we just have to normalize the previous state in order to obtain the excited state |G5 + 1, n d — l)o6---w Up to a sign, 
we have 



\GS+l,n d -l)* b ... v 



1 



^n d (25 + n d - 1)(25- 1) 
E { ( II 4)^25(25-1) |G5+1)4 



i=2,i^i 



EC" 1 )^ II dl)(V2S \GS + l)t jXi ~V2S^2 \GS+l)± iXj : 

j=2 1=1,1^3,1 



and 



c\d a \GS,n d )^ b ... v = (A^2)\GS + l,n d -l)^ b ... v 



3. Tables of Clebsch-Gordan coefficients 

The calculation of the excited states involves the use 
of some Clebsch- Gordon coefficients. We have evaluated 
these quantities explicitely for arbitrary 2S, n d and N, 
following the steps outlined in App endix B 1 . We sum- 
marize our results in Tables IV- VII. In the tables, all the 



coefficients are assumed to be under the sign of the square 
root. For instance, — (25 — 1) corresponds to — V25 — 1. 
The states should be divided by the normalization factor 
VN. Tables [Fv|-|v| are necessary when one electron is 
add ed to the effective impurity, ([1] <g> [25 — 1]). Table 
VII corresponds to the addition of n d electrons. 
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TABLE IV: Clebsch-Gordan coefficients for the process [1] <g> 
[25 - 1] -> [25] [25 - I, I]. The label \{a}aab) indicates 
a state of [25], whereas |{a}aa, b) denotes the highest weight 
state in [25-1, 1] 
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-(25-1) 



TABLE VI: Same as Table [v], but for states with quan- 
tum numbers {a) 2S ~ 2 bc. Notice the degeneracy in [25 — 1, 1]. 
We denote corresponding orthogonal states by |{a}ab, c), and 
|{a}ac, b). 
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numbers (a) 2S ~ 2 bb 



mann, Phys.Rev.Lett. 72, 3262 (1994). 

4 H. v. Lohneysen, J.Phys.: Condens.Matter 8, 9689 (1996). 

5 S. Kambe, S. Raymond, L. Regnault, J. Flouquet, P. Lejay, 
and P. Haen, J. Phys. Soc. Jpn 65, 3294 (1996). 

6 O. Trovarelli, C. Geibel, S. Mederle, C. Langhammer, 
F. M. Grosche, P. Gegenwart, M. Lang, G. Sparn, and 
F. Steglich, Phys. Rev. Lett. 85, 626 (2000). 

7 F. Steglich, P. Gegenwart, R. Helfrich, C. Langhammer, 
P. Hellmann, L. Donnevert, C. Geibel, G. S. M. Lang and, 
W. Assmus, G. Stewart, et al., Z. Phys. B 103, 235 (1997). 

8 F. Steglich, B. Buschinger, P. Gegenwart, M. Lohmann, 
R. Helfrich, C. Langhammer, P. Hellmann, L. Donnevert, 
S. Thomas, A. Link, et al., J. Phys.: Condens. Matter 8, 
9909 (1996). 

9 N. Mathur, F. Grosche, S. Julian, I. Walker, D. Freye, 
R. Haselwimmer, and G. Lonzarich, Nature 394, 39 (1998). 

10 C. Seaman, M. Maple, B. Lee, S. Ghamaty, 
M. Torikachvili, J.-S. Kang, L. Liu, J. Allen, and 
D. Cox, Phys. Rev. Lett. 67, 2882 (1991). 

11 J. Hertz, Phys. Rev. B 14, 1165 (1976). 

12 A. Millis, Phys. Rev. B 48, 7183 (1993). 

13 M. Continentino, Phys. Rev. B 47, 11581 (1993). 

14 T. Moriya and T. Takimoto, J. Phys. Soc. Jpn 64, 960 
(1995). 

15 M. Lavagna and C. Pepin, Phys. Rev. B 62, 6450 (2000). 

16 A. Schroder, G. Aeppli, E. Bucher, R. Ramazashvili, and 
P. Coleman, Phys. Rev. Lett. 80, 5623 (1998). 

17 A. Schroder, G. Aeppli, R. Coldea, M. Adams, O. Stock- 
ert, H. v. Lohneysen, E. Bucher, R. Ramazashvili, and 
P. Coleman, Nature 407, 351 (2000). 

18 P. Coleman, Physica B 259-261, 353 (1999). 











in 


abb 


N 


\a)\{a}bb) 


\b)\{a}ab) 


20 










21 










22 


\{a}abb) 


25 


25-2 


2 










23 










24 


\{a}ab, b) 


25 


2 


-(25-2) 












25 



Q. Si, S. Rabello, K. Ingersent, and J. Smith, Nature 413, 
804 (2001). 

D. Grempel and Q. Si, cond-mat/0207493. 
S. Burdin, M. Grilli, and D. Grempel, cond-mat/0206174. 
J. Gan, P. Coleman, and N. Andrei, Phys. Rev. Lett. 68, 
3476 (1992). 

C. Pepin and M. Lavagna, Phys. Rev.B 19, 12180 (1999). 
P. Coleman, C. Pepin, and A. Tsvelik, Phys. Rev. B 62, 
3852 (2000). 

P. Coleman, C. Pepin, and A. M. Tsvelik, Nucl. Phys. B 
586, 641 (2000). 

26 P. Anderson, Phys. Rev. B 164, 352 (1967). 

27 P. Nozieres, J. Phys. (Paris) 37, CI (1976). 

28 A. Hewson, The Kondo problem to heavy fermions (Cam- 
bridge University Press, 1993). 

29 P. Nozieres and A. Blandin, J. Phys. (Paris) 41, 193 (1980). 

30 N. Andrei and C. Destri, Phys. Rev. Lett. 52, 364 (1984). 

31 P. Wiegmann and A. Tsvelik, JETP Lett. 38, 403 (1983). 

32 D. Cox and A. Ruckenstein, Phys. Rev. Lett. 71, 1613 
(1993). 

33 M. Hammermesh, Group Theory and its application to 
Physical problems (Addison Wesley, 1962). 

34 H. Jones, Group Representations and Physics (Institute of 
Physics, 1990). 

35 J. Cornwell, Group Theory in Physics (Academic, London, 
1986). 

36 D. Lichtenberg, Unitary symmetry and elementary parti- 
cles (Academic Press, New York, 1970). 

37 b. Coqblin and J. R. Schrieffer, Phys. Rev. 185, 847 (1969). 

38 S. Okubo, J. Math. Phys. 18, 2382 (1977). 

39 O. Parcollet, A. Georges, G. Kotliar, and A. Sengupta, 
Phys. Rev. B 58, 3794 (1998). 

40 A. Jerez, N. Andrei, and G. Zarand, Phys. Rev. B 58, 3814 
(1998). 

41 I. Affleck and A. W. W. Ludwig, Phys. Rev. B 48, 7297 
(1993). 

42 I. Affleck and A. W. W. Ludwig, Phys. Rev. Lett. 67, 161 
(1991). 

43 I. V. Schensted, A course on the application of Group the- 
ory to Quantum Mechanics (NEO Press, Maine, 1976). 

44 J. Chen, P. Wang, Z. Lii, and X. Wu, Tables of the 
Clebsch-Gordan, Racah, and subduction coefficients of 
SU(n) Groups (World Scientific, Singapore, 1987). 



27 



TABLE VII: Some of the CG coefficients for the product of k = rid electrons, and an effective impurity [25 — 1], ([l fe ]<8> [25— 1] — > 
[25, l fc_1 ] [25 - 1, l fe ] • • • ). Here we only keep the coefficients for the representations [25, l fc_1 ], and [25 - 1, l k ] (last row). 

States from [l fe ] are denoted by a column of labels, |:}; those from [25 — 1] are denoted \a), and are states close to the highest 
weight state, denoted |o). The states in [25, l fe_1 ] are labeled \ab, c,d, . . . ) and those in [25 — 1, 1 ], \a, b, c, • • • ) 
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